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{fu  Akstrect  f  // 

In  the  analysis  of  the®-^ tracking  filter,  it  is  normally  assumed  that  the  tracking 
filter  and  data  source  operate  in  synchronism  at  a  constant  data  rate.  An  analyti¬ 
cal  solution  is  obtained  for  the  case  in  which  the  tracking  filter  and  data  source 
operate  asynchronously,  thus  violating  the  standard  assumptions.  In  order  to  com¬ 
pensate  for  the  asynchronous  operation  of  the  filter,  the  technique  of  time  correc¬ 
tion  is  used  to  adjust  the  measured  data  point  via  the  estimated  velocity  in  order 
to  approximate  the  synchronous  operation  of  the  filter  and  the  data  source.  A  spe¬ 
cific  example  is  given  in  which  the  influence  of  time  correction  on  the  performance 
of  the  altitude  tracker  is  evaluated  in  terms  of  the  filter  performance  for  extended 
time-interval  position  prediction  as  is  used  in  Conflict  Alert.  It  is  shown  that 
errors  in  the  predicted  positions,  on  the  order  of  a  few  thousand  feet,  can  be  intro¬ 
duced  if  the  time-correction  process  is  not  incorporated  into  the  altitude  tracking 
algorithm.  In  addition,  it  is  shown  that  the  magnitude  of  the  errors  increases 
significantly  as  the  altitude  change  rate  increases,  while  for  an  altitude  tracker 
with  time  correction  the  position  prediction  performance  is  essentially  constant 
for  all  altitude  change  rates  of  practical  significance.  The  computational 
resources  required  for  the  implementation  of  the  time-correction  process  are 
inconsequential  as  compared  to  the  total  requirements  of  the  tracking  functions. 
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EXECUTIVE  SUMMARY 


The  basis  for  all  advanced  air  traffic  control  functions  is  the  ability 
to  predict  the  future  position  of  an  aircraft.  An  important  part  of  the 
prediction  process  is  the  calculation  of  the  future  altitude  of  an  aircraft, 
which  is  performed  by  the  altitude  tracking  function.  It  is  the  objective  of 
the  present  study  to  evaluate  the  performance  of  the  altitude  tracking  func¬ 
tion  and,  in  particular,  to  evaluate  the  influence  of  an  assumption  which  has 
been  made  in  the  design  of  the  altitude  tracker;  namely,  that  all  altitude  data 
can  be  assumed  to  have  been  received  at  the  center  of  the  tracking  cycle.  In 
reality,  however,  the  altitude  data  may  have  been  received  up  to  3  seconds 
before  or  after  the  time  at  which  it  was  assumed  to  have  been  received. 

Although  this  assumption  apparently  results  in  some  minor  saving  of  computa¬ 
tional  resources,  it  introduces  an  additional  source  of  error  in  the  computa¬ 
tions  of  the  altitude  tracker.  It  is  shown  in  this  study  that  the  use  of  the 
above  assumption  may  result  in  an  additional  error  of  up  to  a  few  thousand  feet 
in  the  predicted  position,  as  computed  120  seconds  into  the  future,  and  these 
errors  are  found  in  certain  cases  of  practical  significance.  In  particular, 
it  was  found  that  the  magnitude  of  the  errors  in  the  predicted  position 
increases  as  the  altitude  change  rate  increases,  and  that  for  change  rates  of 
practical  significance,  the  error  in  the  predicted  position  may  exceed  the 
vertical  separation  standards. 

Elimination  of  the  errors  in  the  predicted  position  induced  by  the  above 
assumption  was  found  to  be  possible  using  a  technique  known  as  time 
correction,  in  which  the  estimated  velocity  is  used  to  adjust  (or  correct)  the 
measured  position  for  the  difference  between  the  actual  time  of  measurement 
and  the  assumed  time  of  measurement  as  used  by  the  altitude  tracker.  Using 
the  time-correction  process,  it  was  found  possible  to  significantly  reduce  the 
errors  in  the  predicted  altitude.  Also,  the  performance  of  the  altitude 
tracker  would  now  be  essentially  invariant  for  all  altitude  change  rates  of 
practical  significance.  Not  only  would  the  prediction  errors  be  significantly 
less  than  the  vertical  separation  standards,  but  the  elimination  of  the 
variation  in  performance  with  the  change  rate  would  mean  that  special  consid¬ 
eration,  via  additional  programing  logic,  would  no  longer  have  to  be  used  in 
the  case  of  high-transition-rate  aircraft.  Since  the  computational  resources 
required  to  implement  the  time-correction  process  are  minimal,  it  is  recom¬ 
mended  that  this  procedure  be  used  in  the  altitude  tracking  function.  Not 
only  is  the  magnitude  of  the  errors  sufficient  to  warrant  this  change,  but  the 
resulting  increase  in  the  accuracy  of  the  predicted  position  may  very  well,  in 
fact,  result  in  a  net  decrease  in  the  computational  requirements  if  fewer 
track  pairs  need  to  be  subjected  to  the  more  detailed  computations  required  to 
determine  if  a  conflict  exists.  It  is  highly  probable  that  the  computational 
resources  required  for  the  time-correction  process  are  significantly  less  than 
the  other,  less  satisfying  approaches  which  might  be  taken.  Also,  the  possi¬ 
bility  of  a  reduction  in  the  computational  resources  for  the  above  reason  is 
an  added  advantage. 
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1 .  INTRODUCTION 


The  a-0  tracking  filter  is  a  widely  used  technique  for  performing  the  operation 
of  numerical  differentiation  to  obtain  velocity  estimates  from  noisy  position 
measurements.  The  simplicity  of  the  algorithm  and  the  limited  computational 
requirements  have  resulted  in  the  use  of  this  filter  in  many  practical  sit¬ 
uations  and,  as  a  consequence,  extensive  analytical  studies  have  been  made  of 
the  a-0  filter  (e.g.,  references  1  to  16).  In  virtually  all  of  the  studies 
which  have  been  performed  to  date,  it  has  been  assumed  that  the  data  are  ob¬ 
tained  at  a  constant  rate.  In  general,  however,  this  is  an  unrealistic 
assumption  because  even  for  a  surveillance  radar  rotating  at  a  constant  rate, 
the  fact  that  the  targets  are  moving  means  that  the  time  interval  between 
position  measurements  will  not  be  constant.  The  reason  for  this  is  simply  the 
fact  that  a  moving  target  will  not  necessarily  be  at  the  same  angular  location 
with  respect  to  the  antenna  so  that  while  the  average  data  rate  will  stay  con¬ 
stant,  the  actual  time  between  samples  will  vary.  As  a  result,  most  practical 
situations  do  not  meet  the  assumption  of  a  constant  time  interval  between  data 
points  on  which  most  <*-0  filter  analyses  are  based.  One  particular  study  in 
which  this  assumption  was  not  made  is  in  the  work  by  Cantrell  (references  13 
and  14),  and  it  is  the  objective  in  the  present  study  to  show  how  the  results 
obtained  by  Cantrell  can  be  applied  to  the  analysis  of  the  altitude  tracking 
function  required  for  Conflict  Alert  (references  17  and  18). 

For  the  purposes  of  en  route  air  traffic  control  there  is  an  additional 
reason,  beyond  that  arising  as  a  result  of  moving  targets,  why  the  data 
samples  will  not  be  synchronized  with  the  operation  of  the  tracking  algorithm. 
Since  a  particular  air  traffic  control  center  may  have  from  10  to  15  different 
sensors  providing  surveillance  information,  it  is  obvious  that  the  tracking 
algorithm  could  not  operate  synchronously  with  all  at  the  same  time.  As  a 
result,  the  tracking  algorithm  is  scheduled  to  operate  at  fixed  time  inter¬ 
vals  and  is  to  process  the  surveillance  data  which  have  been  received  since 
the  previous  operation  of  the  tracking  algorithm  (reference  17).  The  specific 
purpose  of  this  study  is  to  demonstrate  the  consequences  of  assuming  that  the 
position  measurements,  in  this  particular  case  the  altitude  data,  and  the 
tracking  filter  operate  in  a  synchronous  manner  when,  in  fact,  this  is  not 
true.  If  the  situation  above  is  recognized,  then  it  is  possible  to  compensate 
for  the  asynchronous  operation  of  the  tracking  filter  and  data  source  by  using 
the  estimated  velocity  to  adjust  the  measured  data  to  compensate  for  the  dif¬ 
ference  in  time  between  the  filter  operation  and  the  actual  measurement  time. 

In  using  such  a  procedure,  known  as  time  correction,  the  degree  of  success  is 
dependent  on  the  ability  of  the  tracking  algorithm  to  provide  accurate  velocity 
estimates  and  the  degree  to  which  the  true  target  trajectory  can  be  expressed 
as  a  first-order  function  of  the  time  difference.  An  explicit  quantitative 
analysis  is  given  which  will  allow  a  comparative  study  to  be  made  between  a 
tracking  algorithm  in  which  the  time-correction  process  is  used  and  one  in 
which  it  is  not  used.  The  performance  statistics  of  interest  in  this  study 
will  be  the  variance  of  the  velocity  estimates  and  the  accuracy  of  the  extended 
time-interval  position  prediction,  both  of  which  are  of  considerable  importance 
in  determining  the  ability  of  the  tracking  algorithm  to  support  functions 
such  as  Conflict  Alert. 
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2.  MATHEMATICAL  ANALYSIS 


2.1  DEFINITION  OF  THE  a  -  6  TRACKING  FILTER. 

The  a  -  8  Cracking  algorithm  is  a  recursive  procedure  which  performs  the  oper¬ 
ations  of  position  smoothing,  position  prediction,  and  numerical  differen¬ 
tiation  for  velocity  estimation  and  is  specified  by  the  equations, 

Z8(k)  -  Zp(k)  +  a (Zm(k)  -  Zp(k)) 

Zyik)  *  Zv(k)  +  (8/T>(ZmCk)  -  ZD(k) ) 

(1) 

Zp(k+1 )  -  Zs(k)  ♦  TZv(k) 

where  Z8(k)  *  smoothed  position  at  the  k^  time  epoch 

Zy(k)  ■  velocity  estimate 
Zp(k)  “  predicted  position 
Zm(k)  *  measurement  position 

T  =  sampling  period  (assumed  constant) 
a,  8  *=  smoothing  constants. 

For  the  purposes  of  the  tracking  algorithm  per  se,  it  is  only  necessary  to 
predict  the  future  position  of  the  target  one  time  interval  into  the  future; 
however,  for  the  purposes  of  advanced  air  traffic  control  functions,  it  is 
necessary  to  make  position  predictions  much  further  into  the  future  so  that 
an  extended  time-interval  position  prediction  will  be  defined  as 

Zp(k,T' )  -  Z8(k)  +  T'  Zv(k)  (2) 

in  which  the  time  interval  T'  is  arbitrary.  The  accuracy  of  the  extended  time- 
interval  position  prediction  is  dependent  on  the  accuracy  of  the  tracking 
filter  outputs,  Z8  and  Zv,  and  also  on  the  degree  to  which  the  actual  flight- 
path  follows  the  constant  velocity,  straight-line  assumption  inherent  in  (2). 

The  algorithm  as  defined  by  (1)  assumes  the  sequence  of  events  as  illustrated 
in  figure  1  with  all  computations  and  measurements  being  coincident  with  the 
epoch  times.  In  an  asynchronous  multisensor  environment,  however,  data  may  be 
received  at  any  time  between  the  operations  of  the  tracking  algorithm.  In  such 
cases,  it  is  necessary  to  assume  a  reference  time  for  the  smoothing  and  pre¬ 
diction  process  which  may  not  necessarily  be  the  time  of  operation  of  the 
tracking  algorithm  or  the  time  of  receipt  of  the  measurement  datum.  In  the 
case  of  the  en  route  portion  of  the  National  Airspace  System,  the  tracking 
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function  operates  at  a  fixed  rate,  not  necessarily  that  of  the  sensor,  with 
the  computational  time  taken  as  the  midpoint  of  the  tracking  cycle  operation 
(reference  17).  The  operation  of  the  tracking  algorithm  is  illustrated  in 
figure  2.  The  smoothing  and  prediction  process  is  assumed  to  use  the  center 
of  the  tracking  cycle  as  the  reference  time,  thus  predicting  from  the  center 
of  the  present  cycle  to  the  center  of  the  succeeding  cycle.  As  illustrated  in 
figure  2,  measurement  data  may  not  be  received  at  the  reference  time  used  by 
the  tracking  algorithm,  and  if  this  is  the  case,  then  the  estimated  velocity 
from  the  previous  cycle  may  be  used  to  move  the  data  point,  either  forward 
or  backward  in  time,  to  make  it  appear  as  though  the  measurement  datum  was 
received  in  synchronism  at  the  center  of  the  cycle.  This  process  is  known 
as  time  correction. 


z8(k-l) 

Z8(k) 

Z8(k+1) 

Zv(k-1 ) 

Zv(k) 

Zv(k+1) 

Zp(k-l) 

Zp(k) 

Zp(k+1 ) 

Zm(k-1) 

Zm(k) 

Zm(k+1) 

—| - 1 - h- 

(k-l)T  kT  (k+1 )T 


FIGURE  1.  TIME  SEQUENCE  OF  EVENTS  FOR  SYNCHRONOUS  TRACKING  FILTER  OPERATION 
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FIGURE  2.  TIME  SEQUENCE  OF  EVENTS  FOR  ASYNCHRONOUS  TRACKING  FILTER  OPERATION 
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In  this 

case  the 

ZB(k) 

Zy(k) 

where 

AT(k) 

smoothing  equations  are 

-  Zp(k)  +  a  (^(k)  ♦  AT(k)  Zv(k-1 )  -  Zp(k)) 

-  Zv(k-1)  +  O/T)  (Z^k)  +  AT(k)  Zv(k-1)  -  Zp(k)), 


-  kT  -  Tn(k), 


(3) 


(4) 


with  Tm(k)  being  the  actual  time  at  which  the  position  measurement  was  made. 


As  a  result  of  the  time-correction  process  it  is  not  even  necessary  for  data 
to  be  received  every  cycle,  since  if  no  datum  is  received  in  a  particular 
cycle,  the  track  (or  assumed  trajectory)  is  simply  predicted  ahead  to  the 
center  of  the  next  tracking  cycle.  (The  opposite  case  in  which  multiple 
measurements  are  received  within  one  cycle  wi!l  not  be  considered,  since  it 
does  not  correspond  to  reality  and  additional  consideration  would  have  to  be 
given  to  the  exact  usage  being  made  of  the  multiple  measurements.) 


Via  the  process  of  time  correction  just  described,  it  has  been  shown  how  it  is 
possible  for  the  tracking  algorithm  to  operate  at  a  fixed  cyclic  rate  and  yet 
the  measurement  data  which  are  used  by  that  algorithm  may  be  obtained  in  a 
totally  asynchronous  manner,  even  to  the  extent  of  being  obtained  at  a  different 
data  rate.  The  multiple  sensor  environment  of  the  en  route  air  traffic  con¬ 
trol  system  meets  the  conditions  just  described.  It  should  be  noted  that  if 
measurements  are  obtained  asynchronously  and  the  time-correction  process  is 
not  used,  then  this  is  equivalent  to  the  introduction  of  an  error  equal  to 
the  difference  between  the  measured  position  and  the  true  position  at  the  time 
the  measurement  should  have  been  made  if  the  requirement  for  synchronism 
between  the  data  source  and  the  tracking  algorithm  had  been  fulfilled.  The 
elimination  or  omission  of  the  time-correction  process  will  introduce  an 
additional  source  of  error  into  the  tracking  algorithm  which  is  unnecessary  if 
the  time  of  receipt  of  the  measured  position  is  known. 


2.2  VARIANCE  REDUCTION  RATIOS  FOR  q -  P TRACKING  FILTERS  USING  TIME  CORRECTION 

The  performance  of  the  a-p  tracking  filter  is  usually  expressed  in  terms  of 
the  variance  reduction  ratios  which  are  the  ratios  of  the  error  variances  at 
the  output  of  the  filter  to  the  variance  of  the  errors  at  the  input  of  the 
filter.  The  variance  reduction  ratios  describe  the  performance  of  the  tracking 
filter  in  a  steady-state  situation  in  which  all  transients  have  decayed.  If 
transient  errors  are  present,  such  as  at  the  start  of  a  maneuver,  then  errors 
significantly  larger  than  those  discussed  in  this  paper  will  be  present.  It 
can  be  shown,  however,  that  the  transient  error  for  constant  velocity  targets 
will  eventually  decay  to  zero  for  the  tracking  filter  regardless  of  whether  or 
not  time  correction  is  used.  Computation  of  the  variance  reduction  ratios 
will  be  facilitated  if  the  tracking  algorithm  equations  are  expressed  in  the 
matrix  form: 
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Z8(k) 

1-a 

T(  1+aAT/T-a) 

Z8(k-1) 

a 

- 

+ 

Zv(k) 

-B/T 

(1+8AT/T-8) 

Zv(k-1) 

B/T 

or, 


(5) 

2*00 


Z(k)  -  A(T,  AT)  Z(k-l)  +  B(T)  (u(k)  +  w(k)) 


(6) 


where, 


Z(k) 


Z  s  (  k ) 
Zv(k) 


A(T,  AT) 


1-  a  T( 1+aAT/T- a) 
-6/T  (1+8AT/T-6) 


B(T) 


a 

B/T 


and  the  measurement  datum,  Zm(k),  is  expressed  as  the  sum  of  a  true  deter-  2 
ministic  component,  u(k),  and  a  random  error  component,  w(k),  with  variance  °w 
which  will  be  assumed  to  be  white  stationary  noise  representing  the  measure¬ 
ment  error. 


The  noise  response  of  the  filter  is  obtained  in  terms  of  the  covariance 
matrix  for  the  errors  at  the  filter  output,  and  this  response  is  given  by 
(reference  19) 

2 

P(k+1 )  *  A(T,AT)P(k)A'(T,AT)  +  B(T)  %  B’(T),  (7) 

2 

where  ow  is  the  variance  of  the  input  noise.  All  of  the  coefficients  in  (7) 
are  constant  with  the  exception  of  AT  which  is  the  random  time-correction 
factor.  Cantrell  has  shown  that  in  the  case  where  matrices  A  and  B  are  random 
variables  which  are  identically  distributed  and  independent  from  sample  to 
sample,  that  the  covariance  matrix  is  given  by  (reference  13) 


2 

P(k+1 )  -  A(T,  AT)P(k)A'(T,  AT)  +  B(T)°WB'(T)  (8) 

where  the  bar  denotes  the  expected  value  (averaged  over  the  random  variable  of 
interest,  in  this  case  AT). 
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To  solve  for  the  variance  reduction  ratios,  A(T,  AT)  and  B(T)  are  used  in  (8) 
with  the  resulting  equations  then  being  averaged  over  AT.  By  performing  the 
required  operations  and  noting  that  in  the  steady-state  case 

P(k+1 )  -  P(k) ,  (9) 


then  (8)  becomes,  after  some  rearranging,  and  assuming  that  E(AT) 


a (2-a) 

-2T( 1-a) 

2 

-T2(l-2a+a2(l+oJT/T2))‘ 

8(l-a)/T 

28-2a8+a 

-T(l-o-B+aB(l+oJT/T2)) 

-(6/T)2 

28(1-6) /T 

2B-82( l+a2x/T2) 

B 

’a2 

aB/T 

(B/T)2 

2 

a  » 
w 

“  0, 
P 


88 


V8 


1 VV 


(10) 


wh4re  a2T  *  E(AT2) 


Pss  =  steady-state  variance  of  the  smoothed  position,  Zs(k) 
pvs  *  steady-state  covariance  of  Zv(k)  and  Z8(k) 
and  Pvv  =  steady-state  variance  of  Zv(k). 


Solving  these  equations  simultaneously  gives 

Ks  =  =  (2a2-3a8+20)/A 

Kvs  "  pvs/°w  *  8(2a-8)/(TA)  (11) 

Ky  =  Pvv/a2  =  2(B/T)2/A 

with  A  *  a (4-2a-0 )  -  2aJT(0/T)2 

where  Ks,  KyS,  and  Ky  are  the  normalized  variance  reduction  ratios  with  respect 
to  the  input  noise.  In  the  case  where  ofa  m  0,  these  equations  reduce  to  2 
the  results  found  elsewhere  (references  3,  12-16,  20).  Since  the  factor  a^T 
tends  to  reduce  the  value  of  the  denominator  in  the  variance  reduction  ratios, 
it  would  appear  that  the  time-correction  factor  would  actually  result  in  an 
increase  in  the  noise  at  the  output  of  a  tracking  filter  in  which  time  cor¬ 
rection  is  used,  but  as  will  be  shown  in  the  following  section,  this  is  not 
the  case. 


In  the  case  of  the  predicted  position,  Zp(k,T')  given  by  (2),  the  variance 
reduction  ratio  can  be  expressed  in  terms  of  the  variance  and  covariance 
reduction  ratios  as 


(12) 


Kp(T')  -  Kfl  +  21'Kvs  +  (T’)2Kv 
or  equivalently, 

Kp(T')  -  (2a2-3aB+28  +  2(T'/T)B(2a-B)  +  2(BT7T)2)/a 


(13) 


where  it  is  explicitly  denoted  that  the  variance  reduction  ratio  for  the 
predicted  position  is  a  function  of  the  time  interval  over  which  the  pre¬ 
diction  is  made.  In  the  case  of  a  single-cycle  prediction  (13),  reduces  to 


Kp  -  (2a2  +  a 8  +  2  8)/i 


(14) 


which  is  the  result  usually  found  in  the  literature  (when  a *  0). 


2.3  APPLICATION  OF  THE  ASYNCHRONOUS  VARIANCE  REDUCTION  RATIOS  TO  ALTITUDE 
TRACKING . 

The  equations  derived  above  can  be  applied  to  the  analysis  of  an  a  - B  tracking 
filter  as  used  for  the  en  route  altitude  tracking  function  (reference  17). 

For  this  purpose,  the  time-correction  process  is  not  used  in  order  to  minimize 
the  computational  requirements,  and  it  is  implicitly  assumed  that  the  error 
introduced  by  neglecting  time  correction  is  justifiable  in  terms  of  simplifica¬ 
tion  of  the  algorithm.  If  the  time-correction  procedure  is  not  performed  in 
an  asynchronous  situation,  then  this  is  equivalent  to  the  introduction  of  an 
error  equal  to  the  difference  between  the  measured  position  and  the  true 
position  at  the  time  at  which  the  filter  assumes  the  measurement  to  have  been 
made.  If  the  target  is  moving  at  a  constant  true  velocity  Zy,  then  the  error 
which  is  introduced  is  equal  to  ATZy  so  that  the  errors  at  the  input  to  the 
filter  can  be  considered  as  two  additive  errors  as  illustrated  in  figure  3. 

The  error  dZ  will  be  assumed  to  arise  as  a  consequence  of  the  quantization  of 
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FIGURE  3.  ILLUSTRATION  OF  INPUT  DATA  ERRORS  FOR  ASYNCHRONOUS  FILTER 
WITHOUT  TIME  CORRECTION 
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Che  Crue  position,  Z?.  It  will  also  be  assumed  that  these  errors  are  white 
and  stochastically  independent.  The  case  in  which  time  correction  is  used  is 
illustrated  conceptually  in  figure  4.  As  seen  in  this  figure,  the  time- 
correction  process  is  a  feedback  loop  in  which  the  estimated  velocity  is 
multiplied  by  AT  to  form  a  corrected  input.  A  second  noise  source  is  also 
needed  in  this  case  to  account  for  the  fact  that  time  is  also  quantized 
so  that  instead  of  the  error  being  ATZy,  it  is  now  ATqZy  where  ATq  is  the 
time-quantization  unit. 

The  performance  of  the  tracking  filter,  in  the  case  where  the  only  errors 
are  those  discussed  above,  can  be  written  in  terms  of  the  appropriate  variances 
and  variance  reduction  ratios.  For  example,  in  the  case  of  the  variance 
of  the  velocity  errors,  the  filter  performance  without  time  correction  is 


where 


Pvv-Kv(°A2Z  +  ZS  %T>. 

Ky  =  2(8/T)2/a(4-2a-8) 

2 

°AZ<e  variance  measurement  quantization  errors. 


0  AT  *  var*-ance  of  AT, 

and'  in  the  case  in  which  time  correction  is  used 


a-/8  FILTER 


A  Ta  Z 
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FIGURE  4.  ILLUSTRATION  OF  INPUT  DATA  ERRORS  FOR  ASYNCHRONOUS  FILTER 
WITH  TIME  CORRECTION 


1 


rvv 


where  Kv  is  given  by  (11)  and 
2 

°AT  *  variance  of  time-quantization  errors. 

q 


(16) 


Similar  results  can  be  specified  for  the  variances  of  the  smoothed  and  pre¬ 
dicted  positions  and  for  the  covariance  between  the  smoothed  position  and  the 
velocity.  If  it  is  assumed  that  all  the  measurement  errors  and  the  time- 
correction  factor,  AT,  are  uniformly  distributed  with  a  mean  of  zero,  then  the 
error  variances  are 

alz~  Az2/12 

°at  "  at2/12  (17) 

and  -  AT2/12, 

q  q 

where  now  AZ,  AT,  and  ATq  are  the  widths  of  the  intervals  in  which  these  quan¬ 
tities  are  contained.  To  evaluate  the  impact  of  time  correction  on  the  perfor¬ 
mance  of  the  tracking  algorithm,  specific  numerical  results  must  be  obtained, 
and  these  results  are  discussed  in  the  following  section. 


3.  NUMERICAL  RESULTS 


The  theory  developed  in  the  previous  section  will  now  be  applied  to  the  specific 
case  of  the  altitude  tracking  function  of  the  en  route  program  (reference  17). 
Since  this  tracking  function  is  presently  implemented  without  time  correction, 
it  would  be  useful  to  determine  if  the  level  of  improvement  which  could  be  ob¬ 
tained  by  the  addition  of  time  correction  would  justify  the  computational  re¬ 
sources  required  to  perform  the  function. 


3.1  VELOCITY  VARIANCE  RATIO. 

To  illustrate  the  performance  differences  between  a  tracking  filter  with  time 
correction  and  one  without,  the  variances  of  the  velocity  errors  will  be  cal¬ 
culated  in  each  case.  Altitude  measurements  are  presently  reported  with  a 
quantization  interval  of  100  feet,  and  it  will  be  assumed  that  the  quantization 
error  is  the  only  source  of  error  in  the  actual  position  measurements.  The 
en  route  tracking  function  operates  on  a  cyclical  basis,  with  a  period  of  6 
seconds  as  the  basic  cycle  time  so  that  AT  =  6s.  (For  the  purposes  of  this 
paper,  the  distinction  between  cycles  and  subcycles  is  immaterial.)  In  the 
case  of  the  Air  Traffic  Control  Radar  Beacon  System  (reference  21)  most  en 
route  sensors  operate  with  a  basic  scan  period  of  10  seconds,  and  time  is 
measured  with  a  quantization  interval  of  1/2  second.  Using  the  results  from 
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the  previous  section,  the  rstio  of  the  velocity  error  variance  in  the  case 
where  time  correction  is  not  used  to  the  variance  in  the  case  when  it  is  used 
is  given  by 


AZ2«-AT2Z*  1  (18) 

r  »  - =-= —  - 

AZ2*AT  Z„  1  -  ( 0 AT/T)2/(6«(4-2a-0)) 

q  V 

and  for  the  parameter  values  presently  used  (a-0.594,  6*0.25),  the  variance 
ratio,  r,  is  plotted  in  figure  5  as  a  function  of  the  altitude  change  rate,  Zy. 

The  parameter  values  describing  the  present  system  are  such  that  the  ratio  of 
the  variance  reduction  ratios,  which  is  given  by  the  second  factor  in  (18),  is 
only  1.0025,  thus  indicating  that  the  errors  introduced  by  the  use  of  time 
correction  are  insignificant.  However,  as  is  easily  seen  from  figure  5,  the 
errors  in  the  case  where  time  correction  is  not  used  increase  significantly 
as  the  altitude  change  rate  increases.  Thus,  the  increase  in  the  velocity 
errors  at  the  output  of  the  tracking  filter  due  to  the  use  of  the  time-correction 
feedback  loop  are  of  no  consequence  when  compared  to  the  effect  of  the  reduction 
of  the  errors  in  the  input  data  due  to  the  use  of  time  correction.  Examination 
of  tfye  results  for  other  smoothing  parameters  of  practical  significance  showed 
a  very  similar  trend.  It  should  be  noted  that  nominal  altitude  change  rates  on 
the  order  of  2,000  ft/min  are  used  by  large  commercial  transports  (reference  22) 
while  the  maximum  change  rates  for  turbojets  are  on  the  order  of  7,000  to  10,000 
ft/min  (reference  28)  so  that  potentially  significant  results  are  likely  to  be 
obtained  for  situations  of  practical  interest.  Needless  to  say,  this  is  obviously 
the  case  for  military  traffic  which  may  use  significantly  greater  change  rates. 


3.2  ERRORS  IN  THE  PREDICTED  ALTITUDE. 

Advanced  air  traffic  control  functions  require  the  prediction  of  the  future 
position  of  an  aircraft  sufficiently  far  into  the  future  to  allow  for  inter¬ 
vention  in  situations  in  which  this  is  warranted.  Since  the  smoothed  velocity 
forms  the  basis  for  the  position  prediction,  it  is  obvious  that  the  velocity 
errors  should  be  directly  related  to  a  more  operationally  identifiable  measure 
of  significance,  and  the  errors  in  the  predicted  position,  as  given  by  (13), 
were  chosen  for  this  purpose.  For  operational  use,  position  predictions  are 
made  for  time  periods  of  120  seconds  and  150  seconds  with  the  longer  predic¬ 
tions  being  used  as  a  filtering  algorithm  to  select  track  pairs  for  a  more 
complete  analysis  in  which  the  120-second  position  predictions  are  used 
(references  17  and  18).  Since  the  prediction  performance  of  the  tracking 
filter  is  highly  dependent  on  the  smoothing  parameters,  results  were  obtained 
for  both  time  intervals  and  smoothing  parameters  currently  used  (reference  17) 
and  for  an  alternative  set  which  has  been  proposed  (reference  23). 

In  order  to  convert  from  the  variance  of  the  predicted  position  to  a  specific 
level  of  error,  some  probabilistic  criterion  must  be  used  at  which  the  error 
is  to  be  specified.  The  probability  chosen  for  this  purpose  was  the  1-percent 
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level  (with  0.5  percent  in  each  tail).  This  particular  level  was  chosen  as 
being  sufficiently  far  out  on  the  tails  of  the  distribution  to  represent  a 
relatively  large  error  but  not  so  far  out  as  to  represent  an  unreasonable 
error.  Since  several  hundred  aircraft  may  be  tracked  at  the  same  time,  a 
portion  of  which  may  be  involved  in  a  vertical  transition,  the  choice  of  a 
1-percent  probability  level  is  not  at  all  unreasonable  because  under  these 
conditions  it  is  highly  likely  that  the  error  in  some  of  the  predicted 
positions  will  equal  or  exceed  those  calculated  at  the  1-percent  level. 
Calculation  of  the  error  level  from  the  variance  requires  the  assumption  of  a 
specific  probability  distribution  for  the  errors  at  the  output  of  the  tracking 
filter.  It  will  be  assumed  that  the  output  errors  are  Gaussian  distributed 
which  is  a  reasonable  assumption  for  practical  purposes  and  is  considered  to 
be  justified  because  of  the  additive  contribution  of  the  input  errors  and  the 
recursive  additive  equations  defining  the  operation  of  the  filter.  It  must 
also  be  noted  that  the  errors  which  will  be  calculated  can  be  either  positive 
or  negative  with  respect  to  the  predicted  position,  so  that  if  a  region  of 
uncertainty  is  to  be  determined,  then  the  size  of  this  region  is  actually 
twice  the  error  as  calculated.  The  relationship  between  the  calculated  error 
and  the  region  of  uncertainty  is  illustrated  in  figure  6. 

Using  the  assumptions  above,  the  error  in  the  predicted  altitude  of  the  1- 
percent  level  as  a  function  of  the  prediction  interval,  T',  is  given  by 

Ze(T')  -  2.576  ^Kp(T')(oJz+ZJ  oJTq) 
in  the  case  where  time  correction  is  used  and 

Ze(T)  -  2.576  y  K'  <T'>  (oj^zj  oJT  ) 

where  K'  <T' )  -  (2q2-3qp+20  +  2(T '/T)B(2a-p)  +  2(pT'/T)2)/a(4-2q-|»)  (21) 

P 


(19) 

(20) 


in  the  case  where  time  correction  is  not  used.  The  errors  in  the  predicted 
altitude  are  plotted  in  figures  7  and  8  as  a  function  of  the  altitude  change 
rate  for  both  prediction  times  and  for  the  parameter  values  of  interest.  As 
expected,  the  same  trend  as  noted  in  figure  5  is  present  in  these  results, 
namely,  that  the  altitude  prediction  errors  increase  significantly  at  the 
higher  change  rates  in  the  case  in  which  time  correction  is  not  applied  to 
the  tracking  filter.  As  the  results  show,  it  is  quite  possible  for  the 
predicted  altitude  to  be  in  error  by  a  few  thousand  feet  for  altitude  change 
rates  used  by  commercial  aircraft. 

When  the  number  of  assumptions  required  to  obtain  these  results  is  considered, 
it  is  very  likely  that  even  larger  errors  will  be  observed  on  occasion,  and 
this  is  especially  true  of  high-performance  aircraft  (e.g.,  military).  When 
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FIGURE  6.  ILLUSTRATION  ON  RELATIONSHIP  BETWEEN  CALCULATED 
ERROR  AND  REGION  OF  UNCERTAINTY 
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FIGURE  7.  ERROR  IN  THE  PREDICTED  ALTITUDE  FOR  SMOOTHING 
PARAMETERS  OF  a-0.59  AND  6-0.25 


it  is  also  considered  that  the  altitude  errors  in  figures  7  and  8  can  be 
either  positive  or  negative  so  that  the  results  must  be  doubled  to  define  a 
region  which  will  exclude  only  1  percent  of  the  predicted  positions,  the  sig¬ 
nificance  of  the  time-correction  process  becomes  readily  apparent. 


3.3  ESTIMATED  TIME  TO  REACH  A  GIVEN  ALTITUDE. 

Another  computation  which  is  sometimes  required  is  the  time  to  reach  a  given 
altitude.  Since  both  the  smoothed  position  and  estimated  velocity  are  used 
in  this  computation,  it  is  of  interest  to  know  how  the  time-correction  process 
will  affect  the  accuracy  of  the  computation.  If  Za  is  a  specified  altitude, 
then  the  estimated  time  to  reach  this  altitude  is 

t  -  (Za  -  Z8)/Zv  (22) 

where  Zs  is  the  current  smoothed  position  and  Zv  is  the  current  estimated 
velocity.  Using  a  first-order  series  expansion,  the  error  in  t,  denoted  as  6t, 
is  given  by 

6t  -  -  AZS/ZV  -  AZv(Za  -  Zs)/zJ  (23) 


where  AZS  andAZv  are  the  errors  in  the  smoothed  position  and  velocity,  respec¬ 
tively.  Since  it  is  assumed  that  all  transients  have  decayed,  the  output  of  the 
tracking  filter  is  unbiased  so  that  E(AZS)  “  E(AZV)  “  0.  Recognizing  that  the 
errors  in  the  smoothed  position  and  the  velocity  are  correlated,  the  variance  of 
the  time  errors  is  given  by 

Ofit  «o^Zv-2(ks  +  2Kvs(Za  -  Z8)/Zy  +  Kv(Za  -  Zs)2/Z2)  (24) 

where  the  variance  reduction  ratios,  Ks,  Kvs,  and  Ky,  as  well  as  the  input  noise 
variance,  are  chosen  to  fit  the  case  under  consideration.  In  the  case  in 
which  time  correction  is  not  used,  (24)  becomes 

°  ° lz  +  Zy  ° AT)ZV~2  |2a2-3aB+2S  (25) 

+  2B(2a-B)(Za-Zs)/TZv  +  2( (Za-Zs)B/TZv)  j /  a(4-2a-fl) 


and  similarly  for  the  case  in  which  time  correction  is  used. 


The  numerical  results  obtained  using  (24)  are  given  in  figure  9  for  an  altitude 
difference  of  4,000  ft.  It  would  be  expected  from  the  nature  of  the  equations 
that  the  timing  errors  would  be  inversely  proportional  to  the  altitude  change 
rate,  and  this  is  the  case  for  low  change  rates;  however,  as  the  altitude 
change  rate  increases,  a  point  will  eventually  be  reached  at  which  the  accuracy 
of  the  time  computation  remains  constant.  The  reason  for  this  can  be  seen 
in  (25)  and  corresponds  to  the  point  at  which  the  first  two  factors  of  Zy 
cancel  and  the  significance  of  the  remaining  factors  has  been  eliminated.  At 
the  lower  change  rates,  the  difference  between  the  performance  with  and  without 
time  correction  is  insignificant,  but  as  the  change  rate  increases,  the  timing 
errors  with  time  correction  become  totally  insignificant,  while  those  in  the 
case  without  time  correction  remain  at  a  level  of  marginal  significance.  In 
the  latter  case,  the  errors  are  sufficiently  large  in  some  situations  to  cause 
an  incorrect  determination  of  the  proper  tracking  cycle  in  which  the  desired 
altitude  would  be  reached.  For  larger  altitude  differences,  the  timing  errors 
observed  without  time  correction  may  correspond  to  two  or  more  tracking  cycles. 


4.  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


An  analysis  was  performed  to  evaluate  the  influence  of  the  assumption  that  the 
altitude  data  are  always  received  at  the  center  of  the  tracking  cycle  when, 
in  fact,  they  are  received  throughout  the  tracking  cycle.  Although  this 
assumption  does  reduce  the  computational  requirements  of  the  altitude  tracking 
function,  the  need  to  emphasize  improved  performance  may  outweigh  the  rather 
modest  savings  in  computational  resources  achieved  with  this  assumption. 

The  technique  used  for  the  analysis  of  the  influence  of  the  above  assumption 
was  based  on  a  similar  situation  analyzed  by  Cantrell  (reference  13)  in  which 
the  influence  of  a  random  update  time  on  the  performance  of  an  a  - B  filter  was 
examined.  In  the  particular  case  of  interest  in  the  present  study,  the 
effect  of  the  assumption  that  all  altitude  data  are  received  at  the  center 
of  a  tracking  cycle  was  analyzed  by  the  introduction  of  an  additional  source 
of  error  in  the  form  of  time- jitter  at  the  input  of  the  tracking  filter. 

After  it  was  found  that  the  impact  on  performance  could  be  significant, 
especially  for  large  altitude  change  rates,  a  technique  was  examined  for  the 
elimination  of  these  errors.  In  this  technique,  which  is  known  as  time 
correction,  the  estimated  velocity  of  the  target  is  used  to  adjust  (or  compen¬ 
sate)  the  measured  position  for  the  difference  between  the  actual  time  of 
receipt  of  the  datum  and  the  assumed  time  of  receipt  used  by  the  tracking 
filter.  Although  it  is  necessary  to  assume  certain  characteristics  of  the 
target  trajectory  to  use  the  technique,  it  has  been  found  to  work  rather  well 
for  the  elimination  of  timing  errors. 

The  numerical  results  obtained  in  this  study,  using  parameter  values  applica¬ 
ble  to  the  en  route  altitude  tracker,  indicate  that  significant  errors  in  the 
predicted  altitude  are  introduced  when  the  time-correction  process  is  not 
used.  In  situations  of  practical  interest,  the  errors  in  the  predicted  posi¬ 
tion  can  be  on  the  order  of  a  few  thousand  feet  and  occur  at  vertical  veloc¬ 
ities  of  practical  significance  with  the  largest  errors  occurring  at  the 
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FIGURE  9.  TIMING  ERRORS  AT  THE  1-PERCENT  LEVEL 
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highest  transition  rates.  The  results  obtained  in  this  study  should  actually 
be  considered  as  conservative  lower  bounds  to  the  errors  which  -..'e  observed 
in  practice,  because  other  errors  which  also  influence  the  performance  of  the 
tracking  algorithm  have  not  been  included.  Such  errors  include  finite  preci¬ 
sion  arithmetic,  deviations  from  a  linear  trajectory  such  as  would  be  caused 
by  random  wind  gusts  which  could  introduce  an  error  in  the  measured  altitude 
greater  than  one  quantization  interval,  altimeter  instrument  errors 
(reference  27),  random  delays  in  transmission  of  the  data  causing  errors  in 
excess  of  0.5-second  time  quantization,  and  deviations  from  the  assumptions 
used  to  obtain  the  results.  An  example  of  this  latter  case  is  the  fact  that 
it  was  assumed  that  the  time  differences  between  the  time  of  receipt  of  the 
altitude  data  and  the  center  of  the  tracking  cycle  were  uniformly  distributed 
and  stochastically  independent  from  cycle  to  cycle.  While  this  may  be  true 
considering  very  long  sequences  of  data,  it  may  not  be  true  for  the  relatively 
short  periods  of  time  in  which  vertical  maneuvers  occur  and  in  which  the 
timing  errors  may  vary  periodically  as  the  phase  relationship  between  the 
tracking  algorithm  operation  and  the  sensor  changes.  As  a  consequence,  the 
peak  errors  in  a  particular  case  may  be  significantly  larger  than  those  pre¬ 
dicted  by  the  analysis  in  this  study.  For  this  reason,  the  time-correction 
process  may  have  an  impact  significantly  greater  than  the  results  of  this  study 
indicate,  especially  when  it  is  considered  that  it  is  the  occasional  extreme 
errors  which  cause  the  false  alarms  in  the  surveillance  algorithms  which  use 
the  tracking  data.  It  is  unreasonable  to  continue  to  use  an  algorithm  designed 
under  an  assumption  that  leads  to  errors  which  invalidate  the  results  for  the 
intended  purposes  under  which  the  algorithm  was  originally  designed. 

As  a  consequence,  it  is  concluded  that  time  correction  should  be  applied  to 
the  vertical  tracking  algorithm  just  as  it  is  now  applied  to  the  horizontal 
tracking  algorithm.  In  the  case  when  time  correction  is  used,  the  performance 
of  the  altitude  tracking  algorithm  is  essentially  invariant  with  respect  to 
the  altitude  change  rate  for  all  rates  of  practical  significance.  Previous 
studies  of  the  vertical  tracking  algorithm  have  all  identified  the  need  for 
accurate  vertical  velocity  information  (reference  23-26),  and  in  particular, 
problems  have  already  been  identified  in  the  case  of  high-transition-rate 
targets  (reference  23).  It  is  reasonable  to  assume  that  the  time-correction 
process  will  eliminate,  or  at  least  mitigate,  the  problems  associated  with 
high-transition-rate  targets  so  that  this  solution  should  be  incorporated 
before  other  approaches  are  attempted. 

It  should  also  be  noted  that  not  only  will  time  correction  result  in  a  more 
accurate  estimate  of  future  aircraft  position,  but  the  simplicity  of  the 
solution  is  such  that  it  probably  has  the  lowest  computational  requirements 
of  any  possible  solution  for  the  following  reason.  Since  the  majority  of 
en  route  traffic  is  in  level  flight  most  of  the  time  (reference  24),  in  which 
case  special  program  logic  will  make  the  altitude  change  rate  identically  zero 
provided  that  the  altitude  datum  is  in  conformance  with  the  assigned  altitude 
(reference  17),  all  that  is  necessary  to  do  to  determine  if  time  correction 
must  be  done  is  to  check  for  a  nonzero  change  rate,  and  for  the  majority  of 
the  time,  this  is  the  only  computation  which  will  be  required.  It  is  highly 
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likely  that  the  computational  requirements  to  perform  the  time-correction 
process  will  be  so  minimal  that  the  change  in  the  computational  resources  re¬ 
quired  to  perform  the  altitude  tracking  function  will  be  almost  undetectable. 
In  fact,  the  use  of  time  correction  may  cause  a  net  reduction  in  the  total 
computational  requirements  if  the  more  accurate  position  predictions  cause 
fewer  tracks  to  pass  through  the  initial  conflict  detection  filter,  which  in¬ 
cludes  a  coarse  altitude  filter,  thereby  reducing  the  number  of  track  pairs 
which  must  be  subjected  to  the  more  detailed  processing  necessary  to  detect 
a  conflict.  However,  regardless  of  the  computational  requirements,  which  are 
inconsequential  even  assuming  the  possibility  of  no  indirect  benefits  as 
just  discussed,  the  errors  resulting  from  the  lack  of  time  correction  appear 
to  be  such  that  this  change  is  absolutely  essential  to  the  proper  functioning 
of  the  altitude  tracking  algorithm. 

Since  in  the  future  other  advanced  air  traffic  control  functions  may  be 
implemented  which  will  require  altitude  predictions,  it  is  recommended  that 
the  time-correction  process  be  incorporated  into  the  altitude  tracking 
algorithm  in  order  to  insure  that  the  altitude  predictions  are  as  accurate  as 
possible.  Otherwise,  it  may  be  necessary  to  incorporate  less  satisfying  solu¬ 
tions  to  the  problems  caused  by  high-transition-rate  targets  which  may  well 
impose  a  far  greater  computational  burden  than  would  be  imposed  by  the  time- 
correction  process  above  and  would  only  have  the  effect  of  suppressing  the 
altitude  prediction  errors  rather  than  eliminating  the  errors  as  is  done  with 
time  correction. 
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